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Q ' We study the behavior of two diferent models at finite temperature in a D-dimensional space- 

^ I time. The first one is the Xip'^ model and the second one is the Gross-Neveu model. Using the 
^ I one-loop approximation we show that in the X(p'^ model the thermal mass increase with the tem- 
perature while the thermal coupling constant decrese with the temperature. Using this facts we 
D ! establish that in the {Xip'^)D=3 model there is a temperature above which the system can 
develop a first order phase transition, where the origin corresponds to a metastable vacuum. In 
the massless Gross-Neveu model, we demonstrate that for D = 3 the thermal correction to the 
coupling constant is zero. For D ^ 3 our results are inconclusive. 
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1 Introduction 



In the last years, there has been much interest in the nature of the electroweak phase transition. 
The high temperature effective potential in the standard and in the {X(p'^)d=4 models have been 
calculated by many authors, where the contribution from multiloops diagrams has been taking into 
account. Several authors have pointed out the importance to known whether in {X(p^) d=4 model the 
phase transition is of first or second order [|l|. Our interest in these issues was stimulated by some 
results of Ford and Svaiter concerning the thermal dependence of the mass and coupling constant 
in (A</?^)£)=4 model defined in a non-simple connected spacetime ||^. In the aforementioned paper 
these authors studied a neutral scalar field in a D = 4 dimensional spacetime using the one-loop 
effective potential. The cases of trivial and non-trivial topology of the spacelike sections and finite 
temperature were discussed. The temperature and topological dependent renormalized mass and 
coupling constant were derived using the Speer and BoUini, Giambiagi and Domingues analytic 
regularization and a modified minimal subtraction renormalization procedure |Q. In addition 
they have also discussed the possibility of vanishing the renormalized coupling constant in this 
model, as well as the limits of validity of the one-loop approximation. Some calculations studying 
such kind of problems was given recently by Elizalde and Kirsten and also Villareal 0. This last 
author improved the precedent results studying the two-loops corrections to the effective potential 
for scalar fields defined in a spacetime with non-trivial topology of the spacelike sections. 

The two goals of this paper are the following. The first one is to extend the discussion of 
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the massive self-interacting X^p'^ model to an arbitrary D-dimensional spacetime, assuming trivial 
topology of the spacelike sections and to analize temperature effects in a model with asymptotic 
freedom. The second one is to discuss the existence of a first order phase transition in the massive 
(A99^)z)<4 model. 

Besides Yang-Mills theories in D = 4, the other known perturbative renormalizable asymptot- 
ically free theories with fermions are the Nambu-Jona-Lasinio and the Gross- Neveu models 0. In 
the latter, a component fermion field with a quartic self-interaction is assumed. The model is 
perturbatively renormalizable for D = 2 and develops asymptotic freedom. 

Working in a generic D-dimensional spacetime, we first calculate the one-loop corrections to 
the renormalized mass and coupling constant in the Xip'^ model. We obtained that the thermal 
mass increase and the thermal coupling constant decrease with the temperature. Still using the 
one-loop aproximation, the thermal correction to the renormalized coupling constant in the Gross- 
Neveu model is obtained. We demonstrate that in the case D = 3 the thermal correction to the 
coupling constant is zero. For D ^ 3 our results are inconclusive. 

In many papers studying second order phase transition in the Xip'^ model the temperature 
dependence of the coupling constant is neglected. This approach is reasonable since the variation 
of the mass with the temperature is the most important fact for a critical phenomena. In this case, 
it is sufficient to consider the renormalized coupling constant as constant and the thermal mass 
drives the second order phase transition [0 . In this paper we will examine the existence of a first 
order phase transition in the X(f^ model taking into account the thermal dependence of the coupling 
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constant. Note that we will not deal with the system behavior in the neighborhood of a second 
order phase transition since we assume that the tree level mass squared is positive. This fact 
prevents the one-loop approximation to break down at low temperatures since there is no infrared 
divergences associated with vanishing masses. The result of our analysis can be summarized as 
follows: for D < 4, there is a temperature where the effective coupling constant vanishes. For 
temperatures (]~^ > (3~^, the renormalized coupling constant becomes negative and the system 
may suffer a first order phase transition. The effects of the radiative corrections is toward the 
direction of breaking a symmetry. Compare with the electroweak first order phase transition p|. 
We should note that at (3~^ = (3~^ the system is still in an interacting phase. For D < 4, there is 
a temperature where only the effective coupling constant {X{/3) = X — }? f{l3)) vanishes. All the 
higher 2n-points correlation functions do not vanish, therefore the model is not gaussian at the 
temperature f3~^ . This is an important point that was stressed by Weldon 

The study of the dependence of the coupling constant with the temperature in QFT is well 
known in the literature. Many authors have studied such dependence in the A(/9^ model and 
also in a abelian model like QED [|T0|. Instead of using perturbative arguments, the use of the 
renormalization group equations allowed the investigation on the mass and coupling constant 
thermal dependence. Such program was implemented by Fujimoto, Ideura, Nakano and Yoneyama 
rT| . These authors obtained results similar to ours in the A(/9^ model. The behavior of the mass and 



coupling constant with the temperature are opposite, i.e. the renormalized mass increases if the 
system temperature increase as where the coupling constant decreases. If we assume that the one- 
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loop approximation provides trustable results we have the following situation: for temperatures 
above (3~^ the renormalized coupling constant becomes negative. This behavior of the effective 
coupling constant is related to the fact that the model is non-asymptotically free. The growth of 
the renormalized coupling constant at large momenta is translated in our case to the temperature 
growth (in modulus) of this quantity. In D = 3 for temperatures /3~^ > the system develop a 
first order phase transition where the origin is a metastable vacuum. 

In this paper we address only the one-loop approximation. It is not unreasonable to believe 
that our conclusions in the A</?^ model may be limited to this approximation. In fact, the behavior 
of the thermal correction to the coupling constant changes in the two-loops approximation. It was 
been shown by Funakubo and Sakamoto [|l^ that only for low temperatures the behavior of the 
thermal coupling constant remains the same as the obtained in the one-loop approximation. For 
high temperatures » m) the behavior is opposite i.e., the thermal correction is positive. 
Nevertheless this fact does not exclude the possibility of a first order phase transition at low 
temperatures in {X(f^)D=3- A more detailed discussion will appear in a forthcoming paper. 

The paper is organized as follows. In section II we sketch the formalism of the effective 
potential. In section III, the massive self-interacting A(/?^ model is analised. In section IV we 
repeat the calculations in the Gross-Neveu model. Conclusions are given in section V . In this 
paper we use ^ = c = 1. 
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2 The effective action and the effective potential at zero 
temperature. 

In this chapter we will review briefly the basic features of the effective potential associated with 
a real massive self-interacting scalar field at zero temperature. Although the formalism of this 
section may be found in standard texbooks, we recall here its main results for completeness. Let us 
suppose a real massive scalar field (p{x) with the usual Xip'^i^x) self-interaction, defined in a static 
spacetime. Since the manifold is static, there is a global timelike Killing vector field orthogonal to 
the spacelike sections. Due to this fact, energy and thermal equilibrium have a precise meaning. 
For the sake of simplicity, let us suppose that the manifold is flat. In the path integral approach, 
the basic object is the generating functional. 



where I>[(/3] is the functional measure and S[iip] is the classical action associated with the scalar 
field. The quantity Z[J] gives the transition amplitude from the initial vacuum \0,in > to the 
final vacuum \0,out > in the presence of some source J{x), which is zero outside some interval 
[—T, T] and inside this interval was switched adiabatically on and off. Since we are interested in 
the connected part of the time ordered products of the fields, we take the connected generating 
functional as usual. This quantity is defined in terms of the vacuum persistent amplitude 



Z[J] =< O,out\0,in > 




(1) 
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by 

e^wiJ] ^<:0^out\0,in> . (2) 

The connected n-point function Gc{xi,X2, is defined by 

6-W[J] , 

Ge(xi,:r2,..,a:„) - JJ^^^j;jJ^^\j=o- (3) 

Expanding VF[J] in a functional Taylor series, the n-order coefficient of this series will be the 
sum of all connected Feynman diagrams with n external legs, i.e. the connected Green's functions 
defined by eq.(3). Then 

OO 1 „ 

W[J] = E -r / dW.d^Xn G'^\x,,X2....Xn)J{xi)J{x2)..J{Xr,). (4) 
n=0 

The classical field (/Pq is given by the normalized vacuum expectation value of the field 

_ <0,out\(fi{x)\0,in>j 
^"^"'^ ~ 5J(x) ~ <0,out\0,in>j ' ^ ' 

and the effective action r[(/?o] is obtained by performing a functional Legendre transformation 

T[cpo] = W[J] - j d^xJ{x)ifo{.x). (6) 
Using the functional chain rule and the definition of given by eq.(5) we have 

^ = -iM. (7) 

Just as W[J\ generates the connected Green's functions by means of a functional Taylor expan- 
sion, the effective action can be represented as a functional power series around the value tpo = 0, 
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where the coeficients are just the proper n-point functions T^^\xi,X2, ■■,Xn) i.e., 

The coefficients of the above functional expansion are the connected 1 particle irreducible 
diagrams (IPI). Actually, T^^\xi, X2, Xn) is the sum of all IP/ Feynman diagrams with n 
external legs. Writing the effective action in powers of momentum (around the point where all 
external momenta vanish) we have 

r[<^o] = / d'x^-VM + l{d,^yZ[^o] + ..). (9) 



The term V{(po) is called the effective potentially |]T3[ .To express V{(po) in terms of the IPI 



Green's functions, let us write r(")(xi, X2, ■■, x^) in the momentum space: 

r(")(xi,a;2, = -i- / d%d%..d''U2n)%h + k2 + -K) e^^^^^'i+'-'^-^-^f^^Hxi, X2, ..,Xn). 

(10) 

Assuming that the model is translationally invariant, i.e. Lpo is constant over the manifold, we 
have 

„ 00 -r ^ \ 

~ (11) 



r[v.o]= /rf^x^- f(")(o,o,...)(^o)" + .. 

n=l ^ 



If we compare eq.(9) with eq.(ll) we obtain that 

n'/'o) = -EAr^"no,o,..)(¥.o)", (12) 

then is the sum of the all IPI diagrams earring zero external momenta. Assuming that 
the fields are in equilibrium with a thermal reservoir at temperature in the Euclidean time 



formalism, the effective potential V{13, fo) can be identified with the free energy density and can 
be calculated by imposing periodic (antiperiodic) boundary conditions on the bosonic (fermionic) 
fields. 

In the next section using the effective potential we will perform the one-loop renormalization 
of the assuming that the system is in equilibrium with a thermal reservoir at temperature 
Since we are interested to make a paralel with the tricritical phenomena where in the tree 
level approximation with V{(fi) = m?(f^ + X^p'^ + a(f^ predicts the existence of a first order phase 
transition if we allow the coefficient of the quartic term to be negative, we will evaluate the effective 
potential in a very unusual way. Instead of summing the series obtaining a log expression, and 
regularizing the model by introducing an ultraviolet cut-off in the Euclidean momenta, we prefer 
to use the principle of analytic extension in each term of the series. The advantage of this method 
lies in the fact that the dependence of mass and coupling constant with the temperature appear 
in a very straightforward way as well as the paralel with the tricritical phenomena. 

3 The one-loop effective potential in the A(/9^ model at 
zero and finite temperature. 

Let us assume the following Lagrange density associated with a massive neutral scalar field: 

^ = \{d,.^uf - \ml ifl - ^(^^ (13) 
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where ifu{x) is the unrenormahzed field and mo and Aq are the bare mass and bare couphng con- 
stant respectively. We may rewrite the Lagrange density as the usual form where the counterterms 
will appear explicity. Defining the quantities 

ip^ix) = {l + SZ)^^ip{x) (14) 

ml = {m'^ + Sm^){l + 5Z)-^ (15) 
Ao = (A + 5A)(l + 5Z)-2, (16) 
and substituting eq.(14),(15) and (16) in eq.(13) we have 

C = lid.^r - - + \md,^f - \^rn\^ - (17) 

where 5Z, 5m^, and b\ are the wave function, mass and coupling constant counterterms of the 
model. After the Wick rotation, in the one-loop aproximation, the effective potential is given by 
^|: 

V(y;o) = VAV^) + ViA^^) (18) 



where. 



and 



M..)^g^P(^A.g)7^ ^^,^J,^^^,^, (20) 

Before continuing, we would like to discuss one important point. Performing analytic or di- 
mensional regularization, we must introduce a mass parameter /x, in terms of which dimensional 

9 



analysis gives to the field a dimension [(p] = ^1/2(^-2) a^d to the coupling constant a dimension 
[A] = fx'^~^. Mass has dimension of inverse length, i.e. [fj] — [m] — L'^, and the effective potential 
(the energy density per unit volume) has dimension of . 

It is not difficult to extend the results given by eqs.(19) and (20) to finite temperature states. 
After a Wick rotation, the functional integral runs over the fields that satisfy periodic boundary 
conditions in Euclidean time. The effective action can be defined as in the zero temperature 
case by a functional Legendre transformation. Regularization and renormalization procedures 
follow the same steps as in the zero temperature case. Although the counterterms introduced at 
finite temperature are the same as in the zero temperature case, the finite part of the physical 
parameters are temperature dependent. In this situation, since the sign of the thermal correction 
to the coupling constant is negative, the possibility of vanishing the renormalized coupling constant 
appears. 

To study temperature effects we perform as usual the following replacement in the Euclidean 
region: 




(21) 



and 



27rn 



(22) 



where ujn — ^ are the Matsubara frequencies. Defining the dimensionless quantities: 



(23) 
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and 

= a-\ (24) 
the Born terms plus one-loop terms contributing to the effective potential give, 

where, 



and 

Vrr(R (n^\ = , , , , , , 

Is VSttV V /i y y (27r) 



^//(/?,V^o) = -^2^ ^„ Ur^ ~ / T^^A (s.a)- (26) 



The function 



A5^(s, ai, .., Oat) = ^ (ain^ + a2n2 + ... + OArn^ + c^) * (27) 



?ii,n2..n]v=— cxD 



is the inhomogeneous Epstein zeta function ||14||, and finally 



Note that the mass parameter introduced in eqs.(23) and (24) will be used from now on, since 
we must have dimensionless functions when working with analytic extensions. 
Let us define the modified inhomogeneous Epstein zeta function as 

oo 

E^(s,ai,a2, ..Oat) = ^ (airi^ + .. + OArn^ + c^)"^ (28) 

ni,n2vnjv=l 
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Defining the new coupling constant and a new vacuum expectation value of the field (I) (dimen- 
sionless for D = 4), 

^ = (30) 

we rewrite eq.(26) withouth use the definition of the inhomogeneous Epstein zeta function as, 

Vn{PA) = ^''^aY.^-^g^'^ d'k (32) 

s=i ^■^ n=-oo"' [an^ + c^ + k^Y 

To regularize the model we will use a mix between dimensional and zeta function analytic regu- 



larizations. Let us first use dimensional regularization|T5[. Using the well known result, 



{k^ + a?y r(s) ' 2' a 

eq. (32) becomes 



[ d'^k _ 7ri d 1 

J fF + a2)^ Tis) ^' ^ ' 



°° ^-1)^+^ . ,o„ 7r5 d. 



Defining, 



r,,(/3,0) = /i^v^E^^/<^^Wr(^-^) E , (34) 



/(D, s) = fid +l,s) = ^P^^^(s - (35) 



and substituting eqs.(27) and (35) in eq.(34) we obtain. 



Vn{P,^) = fi''y^Y.fiD,s)g'cl^''Af{s - -,a). (36) 
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As we will soon see, the terms s < y are divergent and we will regularize the one- loop effective 
potential using the Principle of the Analytic Extension. Let us assume that each term in the series 
of the one-loop effective potential V{(3, (p) is the analytic extension of these terms, defining in the 
beginning in an open connected set. To render the discussion more general, let us discuss the 
process of the analytic continuation of the modified inhomogeneous Epstein zeta function given 
by eq.(28). For Re{s) > y, the ai, 02, ..oat) converges and represent an analytic function 

of s, so Re{s) > Y is the largest possible domain of the convergences of the series. This means 
that in eq.(36) in the case = 4 only the terms s = 1 and s = 2 are divergent. The term s = 1 
is the divergent one-loop diagram of the connected two-point function and it contributes with a 
quadratic divergence. The s = 2 term is the divergent one-loop diagram of the connected four 
point function, and it contributes to the effective potential with a logarithmic divergence. Using 
a Mellin transform it is possible to find the analytic extension of the modified inhomogeneous 
Epstein zeta function. After some calculations using Kirsten's results [Q, we have: 



(37) 



where K^{z) is the Kelvin function |(T^ 
It is not difficult to show that: 



=1 
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where: 

1 1 f-l)*+i 1 
HD..)^^^L^^y (39) 

If we suppose that D — 4, the model is perturbatively renormahzable and an appropriate choice 
of Sm^ and 5X will render the analj^ic extension of the terms of the series in s in the effective 
potential analytic functions in the neighbourhood of the poles s — 1 and s = 2 respectively. 

The idea to extend the definition of an analytic function to a larger domain (analytic ex- 
tension) and subtract poles was exploited by Speer, BoUini and others. In the method used by 
BoUini, Giambiagi and Domingues, a complex parameter s was introduced as an expoent of the 
denominator of the loop expressions and the integrals are well defined analytic functions of the 
parameters in the region Re{s) > Sq for some Sq. Performing an analytic extension of the expres- 
sion for Re{s) < Sq, poles will appear in the analytic extension and the final expression becomes 
finite after a renormalization procedure. To find the exact form of the counterterms let us use the 
renormalization conditions 

^m'^)l0=o = mV (40) 

and 

^m0)l0=o = A/. (41) 

Since the vacuum expectation value of the field has been chosen to be constant, there is no need 
for wave function renormalization. Substituting eqs.(25),(38) and (39) in eqs.(40) and (41) it is 
possible to find the exact form of the countertems in such a way that they cancel the polar parts of 
the analytic extension of the terms s = 1 and s = 2. Note that we are using a "modified" minimal 

14 



subtraction renormalization scheme where the mass and couphng constant counterterms are poles 
at the physical values of s. It is straighforward to show that both dm'^ and 6X are temperature 
independent. If a model at zero temperature is renormalizable with some counterterms it is also 
renormalizable at finite temperature with the same counterterms. This result was obtained in 
all orders of perturbation theory by Kislinger and Morley [0. In the neighbourhood of the 
poles s = 1 and s = 2, the regular part of the analytic extension of inhomogeneous Epstein zeta 
function has two contributions: one which is temperature independent and that can be absorbed 
by the counterterms and another that is temperature dependent and cannot be absorbed by the 
counterterms. It is clear that the temperature dependent mass is proportional to the regular part 
of the analytic extension of the inhomogeneous Epstein zeta function in the neighborhood of the 
pole s = 1. The same argument can be applied to the renormalized coupling constant. The thermal 
contribution to the renormalized coupling constant is proportional to the analytic extension of the 
inhomogeneous Epstein zeta function in the neighborhood of the pole s = 2. The choice of the 
renormalization point = implies that only the regular part in the neighborhood of the pole 
s = 1 will appear in the renormalized mass. In the next section (where massless self-interacting 
fermion fields are studied) we will show that all the terms of the series of the effective potential 
contribute to the renormalized mass and coupling constant and the sign of the thermal coupling 
constant cannnot be computed for D ^ 3. From the above discussion we can write 

- f(2) {p = 0, /?, A, m) = m^{(3) = + ATu^iP) (42) 
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and 

- f^^\p = 0, A, m) = A(/3) = A + AA(/?), (43) 

where m^(/3) and A(/3) are respectively the temperature dependent renormahzed mass squared and 
couphng constant. It can be directly shown that the thermal contribution to the renormahzed 
mass squared is given by: 

Am\(3) - Am^oo) = —A Y -^KAmnp). (44) 

Using the asymptotic representation of the Bessel function Kn{z) for small arguments 

K„(2) = ir(n)(^)-" ,2^0 n = l,2,.., 

we obtain that at high temperatures the temperature dependent mass squared is proportional to 
A/?~^ [jl9|. The result given by eq.(44) was also obtained by Braden using Schwinger's proper 



time method. The same author also discussed the two-loop effective potential and the problem 
of overlapping divergences where the possibility of temperature dependent counterterms appears. 
Nevertheless these divergences must cancel as it was stressed by Kislinger and Morley . 



Based uppon the same arguments previously used, the thermal contribution to the renormahzed 
coupling constant is given by: 

Q OO 

AA(/?) - AA(oo) = -— A^ Y Ko{mn(3). (45) 

The Bessel function Ko{z) is positive and decreases for ^ > 0. Therefore let us present an 
interesting result: the renormahzed coupling constant attains its maximum at zero temperature 
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= oo) and decreases monotonically as the temperature increases. In other words, the thermal 
contribution to the renormahzed couphng constant AA(/3) — AA(cxd) is negative, and increases in 
modulus with the temperature. The same result was obtained by Fujimoto, Ideura, Nakano and 



Yoneyama using the renormalization group equations at finite temperature |rT|. Once we are 
discussing thermal effects, in the limit of zero temperature the thermal contribution to the mass 
and coupling constant must vanish {T^'^\p = 0,/3 = oo,A,m) = — and (T^^\p = 0,(3 = 
oo,X,m) = —A). This can be easily seen from eqs.(44) and (45). Since the thermal contribution 
to the renormalized coupling constant is negative someone could enquiry: is it possible for the 
renormahzed coupling constant to vanish? Once AA(/3) is O(A^) and we assume D = 4, it is not 
possible to implement such a mechanism for finite temperatures. For D < 4 the renormalized 
coupling constant is not necessarily a small quantity and it can even become a large quantity, due 
to its positive dimension 4 — D in terms of the mass parameter (or using the language of critical 
phenomena, due to its positive dimension 4 — D in terms of the scale ^ where a is the lattice 
spacing). Therefore we conclude that in the neighbourhood of D = 4, the renormalized coupling 
constant A(/3) could vanish only for very high temperatures. As we consider smaller spacetime 
dimensions the temperature where A(/?) vanishes becomes lower and lower. For instance, for D = 3 
we expect to find a finite temperature such that the renormalized coupling constant vanishes. 

We note that there is no discontinuity in the behavior between the cases D = A and D < 4 
as we will see later (see eq.(49)). For D < 4 the model becomes superrenormalizable and only a 
finite number set of graphs need overall counterterms. In the one-loop aproximation for D = 4 
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there are only two divergent graphs and for D < 4 there is only one. This result can be easily 
obtained by investigating eq.(38). In this equation the divergent part of the effective potential is 
given by T(s — y) and for D < A only the s = 1 pole will appear. In other words, for D < A there 
is only finite coupling constant renormalization at the one-loop aproximation. The graph s — 2 
gives a finite and negative contribution to the coupling constant. For > 4 the renormalization 
of the couphng constant is obhgatory (note the presence of the pole in s = 2). Going back to 
the D-dimensional case, the renormalization conditions also are given by eqs.(40) and (41). Using 
the renormalization conditions in eq.(38), we can find the regular part of the analytic extension 
which gives a finite contribution to the renormalized mass squared Am?{D,m, X, /3) and coupling 
constant AX{D, m, A, /?) in a £)-dimensional flat spacetime. We will simplify the notation writing 
Am^{/3) and AA(/3). For even D they are given respectively by: 

and 

where ip{s) = ■^lnT{s). For odd D, the first term between parentesis in eqs.(46) and (47) must be 
replaced by r(l — y)(^)^~^ and r(2 — y)(^)^~^ respectively. The first terms between parentesis 
of eq.(46) and eq.(47) are temperature independent therefore it is possible to isolate the thermal 
contribution to the renormalized mass and coupling constant in a generic D-dimensional spacetime 
in the one-loop aproximation. Using eq.(46) and eq.(47) we obtain the following contribution to 
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the thermal mass and couphng constant respectively: 



,.D-2\ / m \ f-i 



and 



Am - AA(oo) = -!g||(-5-)-^A-„_,(„,,«/3). (49) 

These are among the main results of the paper. Since AA(/?) — AA(oo) < we may have a 



temperature f3~^ where A(/?) vanish for D < A. Our result is different from the Frohlich result |2T 
in which all all the Green's functions of the theory for D > 4 correspond to a free field i.e. the 
model is gaussian at zero temperature above four spacetime dimensions. In our case, the higher 
2n-point functions are not zero as was discussed by Weldon 0. 

Before discussing a existence of a first order phase transition, we would like to point out that 
the investigation of the (Ay?^) 0=4. niodel with a negative bare coupling constant has recently been 
done by Langfeld et al, where an analytic continuation of the model with positive A to negative 
values was presented 0^ . Although several authors claim that the renormalized coupling constant 



of the X(fi model must be positive, a definitive supporting argument is still lacking. Previous 



investigations have been done by many authors |23|. We would like to stress that the sign of the 



renormalized coupling constant is not fixed by the renormalization procedure in the {X(p'^)d=4:- 
Gallavoti and Rivasseau discussed examples with positive bare coupling constant where different 



cutoffs lead to renormalized coupling constants with different signs E 



Going back to the discussion of a first order phase transition, let us define a dimensionless 
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effective potential = -tt, as: 

+ /ligi/i order terms in s. (50) 

In the effective potential all the powers of the field will appear as stated in eq.(38). For 

instance, the term corresponding to the 2s — th power of the field is proportional to 

|(^)*""^f-.(-/3)*='- (51) 

The previous results can be used to demonstrate a first order phase transition in the {\(p^)D=d, 
model. To simplify our discussion let us assume that is possible to truncate the series of the 
effective potential in s = 3. These does not imply the assumption that high order powers of the 
field gives vanishing contributions. They are simply neglected as compared to the leading terms, 
since we are interested in the profile of the effective potential near the origin. The coefficient of 
ip^ is positive (one requires this to ensure that the truncated effective potential is bounded from 
below). For the sake of simplicity, let us also assume that the coefficient of the is constant 
and given by a for both cases = 3 and D = 4. In these cases the leading contributions to the 
effective potential are respectively: 



4(27r)^ ^iV/u2/?n. 
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and 



167r2 :^^\ij?(3n 

+ (i!-^E^o(mn/?))0^ + a0^ (53) 

From the above discussion, for D < 4 we obtain the following profile for the effective potential in 
the neighborhood of the origin. Bellow the temperature the dimensionless effective potential 
has only one global minimum. Heating the system above the temperature (3~^, the renormalized 
coupling constant would become negative and the system can develop a first order phase transition 
since the expectation value of the order parameter changes discontinuously by temperature effects. 
The situation is similar to the Coleman- Weinberg mechanism for massless fields. The effects of 
the quantum corrections is towards the direction of breaking a symmetry. Note the similarity with 
the tricritical phenomena where in the tree level (^(v?) = rn^^"^ + A(/9'^ + CTLp^) the model develop 



a first order phase transition if we allow the coefficient of the quartic term to be negative p5 . 

In a detailed study, using the ring-improved one-loop effective potential, Arnold and Spinosa 
ZB[ showed that even for temperature independent coupling constant, the Ay?'^ model can develop 



at the first sight a first order phase transition. Nevertheless, these authors verified that the con- 
tribution of higher loop corrections dominates over the one-loop ring improved contributions. By 
these reasons, in this approximation they cannot distinguish between a first or a second order 
phase transition. As we discussed in the introduction, the thermal correction to the coupling 
constant if we include high order loops in the effective potential is positive for high temperatures. 
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Nevertheless for low temperatures the effective renormalized coupling constant may become neg- 
ative. In this case we still have a first order phase transition. From the above discussion, we 
have obtained the following result: in the massive Ay?'^ for D < 4 for temperatures above the 
effective potential will develops a local minimum at the origin (a false vacuum) and a global one 
outside the origin. In this case the initial metastable phase may decay to a stable one by nucle- 
ation of bubbles. The temperature is the parameter that drives the first order phase transition. 
Evaluating the ring diagrams Carrington and Takahashi independently obtained in a pure scalar 



model at D = A results which are consonant with ours results in D = 3 |^| 



4 The one-loop effective potential in the massless Gross- 
Neveu model at finite temperature. 



Our purpose throughout this section is to examinate the behavior of the renormalized coupling 
constant in a model involving fermions with a quartic interaction. In two-dimensional spacetime 
{D = 2) the model is renormalizable and ultraviolet asymptotically free. We will consider an 
N-component fermion field where the limit of large N will be investigated. As it was discussed in 
ref.(4), due to the quartic nature of the interaction, it is possible to introduce an ultralocal auxiliar 
scalar field ip which is formally equal to gipip where is the fermionic field, in order to present 
the effective potential of the model. As we did in section II, we suppose that the quantum field is in 
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thermal equilibrium with a reservoir at temperature We will show that for D = 2 and D = A 
inthe one-loop approximation the sign of the thermal correction to the renormalized coupling 
constant cannot be calculated. On the other hand, for D = 3 in the one-loop aproximation the 
thermal correction to the renormalized coupling constant is zero. 
The Lagrange density of the massless model is given by: 

C^i/j, ^p, tp) = ii^-f^df.ip - - gfipij. (54) 

Defining tpo as the vacuum expectation value of tp, i.e. ipo =< 0\(p\0 >=< Olgi/ji/jlO >, the leading 
terms in the effective potential for large N are given by the tree-level graphs plus all one-loop 
graphs, 

H^)4.?-.iv|^(s.o)-/(gj^. (55) 

After a Wick rotation we identify the effective potential as the free energy of the system. At 
zero temperature the model has a spontaneous breakdown of the chiral symmmetry where the 
fermions acquire mass. The symmetry is restored at finite temperature by a second order phase 
transition [^. This result can be obtained by summing the series in the effective potential. Since 
we are interested only in the thermal behavior of the mass and coupling constant instead of repeat 
the well known calculations we will adopt a very unusual road, similar to the previous chapter, by 
regularizing each term of the series in the effective potential before summing up. 

To introduce finite temperature effects we assume that the Grassmannian integration in the 
path integral goes over anti-periodic configurations in Euclidean time. In the effective potential 
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this is equivalent to the replacement given by eq.(21) and 



^-^(^+^)- (56) 



Using eq.(33) and defining f{D, s) by: 

1 (-l)^r(s- 



P(A^) = ^^^^^, (57) 
it is not difficult to show that V{P, tpo) is given by: 

1 OO OO 1 

s=l n=-oo + 2'' 

Note that we are using dimensional regularization in eq. (55) and it is well known that for massless 



fields this technique requires modification in order to deal with infrared divergences [^]. Since we 
are regularizing only a d = D — 1 dimensional integral, this procedure is equivalent to inserting a 
mass into the d dimensional integral. In other words, the Matsubara frequency plays the role of a 
"mass" in the integral, provided we exclude the limit p —* oo, which means that we must restrict 
ourselves to non-zero temperature. 

Again, as in eq.(30), we can define a new field <t> = ^ (no confusion must be done between the 
present auxiliar scalar field and the previous scalar field). Using eq.(24) we obtain 



-1 OO OO -| 

y(/?,0) = -^V + iv^''EMAs)a^-^(#)'^ E 7:rTT^,- m 

The Hurwitz zeta function is defined as 



(r) -I- l )'2s-d ■ 
s=l n=-oo ^ 2> 



OO -| 

C(^,?) = Er^^ (60) 
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for Re{z) > 1 and q ^ 0, —1, .... For q = 1 we recover the usual Riemann zeta function. Defining: 
r(D, .) = piD, s) (C(2. -d^-) + i-ir-'a2s - d, -\) - ^) (61) 
the effective potential can be written as: 

-1 oo 

0) = + iV/.^ s)af-%g<j>r. (62) 

The effective potential is still baddly defined and it will be regularized by the principle of analytic 
extension. The function r{D, s) is valid in the begining in an open connected set of points, i.e. for 
Re{z) > 1. Since we are considering even non perturbative renormalizable models, let us study the 
cases D = 2,3 and 4. We would like to stress that even for the non perturbative renormalizable 
models it is possible to make qualitative predictions and we will regularize and renormalize the 
model in the standard way. A strong argument in favor of the study of the Gross-Neveu model is 
that the non-renormalizability does not appear in the leading approximation for D = 3. 

After the analytic continuation, the effective potential requires a renormalization procedure in 
the points s = 1,2.. The renormalization condition which will fix the form of the counterterm of 
the pole s = 1 is: 

-Q^U=cte = ^^ (63) 

Due to infrared divergences, we must follow Coleman and Weinberg |]l3l and choose the renormal- 
ization point at non- zero 0. In order to evaluate the renormalized effective potential it is necessary 
to use the Hermite formula of the analytic extension for the Hurwitz zeta function given by PD | 

1 0^^^ —z 7/1 

({z,q) = — + ^^ + 2 (g^ + y^)^sin(zarctan^) dy. (64) 

2q^ z — 1 Jo q e'^^^ — 1 
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It is not difficult to show that the thermal contribution to the renormalized coupling constant is, 

oo 

A^(/?) = N^Ji^-^Y.^{D,s){2s){2s-l)g'\^^xr-^, (65) 

s=l 

where it is understood that the polar terms in the summation have been subtracted remaining 
just the regular part of the analytic continuation. The situation is different from the massive Xtp'^ 
model, since we have the contribution of all terms of the series in s and the sign of the thermal 
contribution to the renormalized coupling constant cannot be easily obtained. Nevertheless, for 
sufficiently small g the leading term is 0{g'^). In this case, for D = 3 and using the fact that 
C(0, g) = 1 — g, we obtain that A^f = 0. We found here that there is no thermal correction to the 
coupling constant at least in the one-loop approximation. Note that ^g{l3) is still not well behaved. 
The terms s > ^ are divergent in the low temperature limit (the use of dimensional regularization 
in the begining of the calculations leads to this situation). For s < the high temperature limit 
of the model is problematic due to the well known fact that ultraviolet divergences are worst as 
the spacetime dimension increases. 

5 Conclusions 

In this paper we studied the renormalization program assuming that scalar or fermionic fields are 
in equilibrium with a thermal reservoir at temperature We have attempted to analize the 
consequences of the fact that not only the renormalized mass, but also the renormalized coupling 
constant acquire thermal corrections. 
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It is well known that if we have a one spatial dimension compactified system at a finite tem- 
perature, which has a spontaneous symmetry breaking there are two different ways to restore the 
symmetry. Since the compactification of one spatial dimension gives us the well known mechanism 
of topological generation of mass, it is possible to restore the symmetry by thermal or topological 
effect. There is a very simple way to interpret the origin of the thermal and topological mass and 
couphng constant. The effective potential is not well defined. Using the Principle of the Analytic 
Extension, we regularize the model and the introduction of counterterms remove the principal 
part of the analytic extension, and the model becomes finite. Meanwhile, in the neighbourhood of 
the poles, the regular part of the analytic extension does not vanish. These temperature depen- 
dent regular part around the poles s — 1 and s = 2 (for D — 4) are identified with the thermal 
correction to the mass and coupling constant. 

It was proved that in the \(p^ model, in the one-loop aproximation, the thermal correction 
to the renormalized mass is positive and the thermal correction to the renormalized coupling 
constant is negative. In this case the renormalized coupling constant attains its maximum at zero 
temperature and decreases monotonically as the temperature increases. Since in D = 4, AA(/?) 
is O(A^) it is not possible to vanish the renormalized coupling constant at a finite temperature of 
the thermal bath. For strong couplings {D < 4) there is a finite temperature where this can be 
achieved. For temperatures > (negative coupling constant) the system can develop a first 
order phase transition, where the origin is a false vacuum. 

It is not all clear for us if at D = 4 the system can develop a first order phase transition. We are 
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using the following argument to disregard such possibility. As we discussed in the introduction, in 
the two-loops approximation at high temperatures the thermal correction to the coupling constant 
is positive. The fact that in D = 4 the model has a small zero temperature coupling constant 
eliminate the first order phase transition in D = 4. 

We would like to emphasize that the massive A(/?^ model does not belong to the same univer- 
sality class of the Ising model. It is well known that it is possible to compare the \(p^ model in 
continuous D-dimensional Euclidean space with the Ising model. One lattice formulation can be 
done and the continuum limit of the model (a 0, where a is the lattice spacing) exist if the 
correlation length goes to infinite. This fact implies that at the continuum limit of the lattice 
model the system must suffer a second order phase transition. In other words, close to the critical 
temperature a D-dimensional Ising model has the same correlation functions as those for a field 
theory (A</?^ model) defined in a D-dimensional Euclidean space near the critical temperature. 
Since in the paper we assume that the tree level mass squared is always positive and we found 
that the thermal mass squared is also positive, we are always far from the critical temperature. 
By these reasons the system cannot fall into the universality class of a Ising model. 

The analysis of this paper suggest two possible directions. First, we have to study the decay of 
the metastable ground state in the (A93^)d<4 model evaluating the nucleation rate per unit volume 
in the system. The theory of bubbles nucleation at zero and finite temperature was proposed 
and developed by many authors [0. The basic result is that the probability per unit volume 
per unit time of the metastable vacuum to decay is given by F = Ae~^^'^\ where S{lp) is the 
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Euclidean action of the "bounce" solution which describes tunneling between a metastable and a 
true vacuum. Another possible direction is to examinate if the met ast ability of the system (the 
false ground state) can be eliminated in a more general scalar model. This former subject will be 
presented soon in a forthcomming paper We conclude the paper with some some questions 
which remain to be answered. 

(i) Is the existence of the first order phase transition in {Xip^)d=3 an artifact of our approxi- 
mation? It will be interesting to obtain a non-perturbative argument to demonstrate or disprove 
this fact in a general way. 

(i) Is the series given by eq.(69) Borel summable ? It is well known that the lack of Borel 
summability means that the system is unstable, since the vacuum to vacuum amplitude develop 



and imaginary part p3|. It would be interesting to investigate these questions. 
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